Simple proof of equivalence between adiabatic quantum computation and the circuit 

model 
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We prove the equivalence between adiabatic quantum computation and quantum computation in 
the circuit model. An explicit adiabatic computation procedure is given that generates a ground 
state from which the answer can be extracted. The amount of time needed is evaluated by computing 
the gap. We show that the procedure is computationally efficient. 
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Introduction. — In the effort to realize a quantum com- 
puter (QC), adiabatic quantum computation [H offers 
a promising alternative to the standard "circuit model" 
0, In comparison to the circuit model, AQC allevi- 
ates the need to perform fast quantum logic operations 
and measurements, which is particularly troublesome in 
the context of fault tolerant quantum computation 
In AQC, the answer to a calculation is contained in the 
ground state of a quantum Hamiltonian. By placing a 
system in the ground state of a simple Hamiltonian and 
then adiabatically changing until the desired Hamilto- 
nian is reached, one carries the system into the compu- 
tationally meaningful state. The AQC model was known 
from the outset to be efficiently simulatable by the stan- 
dard model [l|, but for some time researchers won- 
dered whether AQC could efficiently simulate the stan- 
dard model. 

Recently, a number of relatively complex proofs of the 
equivalence between the circuit model and AQC were 
given. One proof showed that AQC using Hamilto- 
nians with long-range five- or three-body interactions, 
or nearest-neighbor two-body interactions with six-state 
particles, can efficiently simulate the circuit model 
This result was soon modified to qubits with two-body 
interactions 0, and then it was shown that AQC us- 
ing qubits with nearest-neighbor two-body interactions 
on a 2D lattice can efficiently simulate the standard cir- 
cuit model [i]. The proofs in Refs. [1, 0, i, H all 
start from a five-body interaction Hamiltonian that arises 
in Kitaev's quantum NP-complete "local Hamiltonians" 
problem [loj . They then require a reduction to two- 
body Hamiltonians, and a proof that the spectral gap 
of the AQC Hamiltonian thus constructed is properly 
lower bounded. Let the number of algorithm steps (num- 
ber of single- and two-qubit gates) be N. The running 
time is 0{N^) with five-body interactions and 0{N^'^) 
with three-body interactions [6|, which was improved to 
0{N^^) with two-body interactions in Ref. [1]. An ad- 
ditional improvement by a factor of N was given in Ref. 

where relatively simple methods were used to pro- 
vide a lower bound on the minimal energy gap. 



Here we provide an alternative, constructive proof of 
the equivalence between the standard circuit model and 
AQC that is physically and mathematically transparent, 
amenable to implementation and yields a running time 
T of order (MN)"^ or better, where M is the number of 
qubits. E.g., in the case of Shor's algorithm for factor- 
ing an L-bit integer using a linear nearest-neighbor qubit 
array [T3|, this translates into T - [(2L + 4)(8L'*)]2 ~ 
256L10 compared to T ~ (SL^^)" ~ W^'^L^'^ using the 
previous 0{N^^) scaling. We do this by setting up an ex- 
plicit Hamiltonian involving at most two-body, nearest- 
neighbor interactions between particles on a 2D lattice. 
Our construction uses the method of ground state quan- 
tum computation (GSQC), which was independently pro- 
posed in Refs. [l3l.[l^ [l5| ar ound the same time as AQC 
and also studied in Ref. Iq . 

In GSQC, one executes an algorithm by producing a 
ground state that spatially encodes the entire temporal 
trajectory of the algorithm, from input to output. This 
requires N times as much hardware but provides some 
robustness against decoherence. GSQC was deliberately 
constructed to simulate the standard model [l^]- How- 
ever, little attention was devoted to the process of reach- 
ing the desired ground state. Here, we marry together 
AQC and GSQC. The result is a formalism supplying an 
explicit Hamiltonian H{s) acting on qubits with at most 
two-body nearest neighbor interactions for any algorithm 
formulated in the circuit model. The initial Hamilto- 
nian H{0) and its ground state are simple. The inter- 
mediate Hamiltonian H{s) (0 < s < 1) has a gap and 
a non-degenerate ground state for all s (the dimension- 
less time). The final Hamiltonian H{1) has a ground 
state containing the solution to the algorithm. Using the 
adiabatic theorem, we provide an upper bound on the 
time needed to reach H{1) while keeping the system in 
its ground state. This bound scales polynomially in the 
number of algorithm steps and qubits; the calculation is 
efficient. 

Single qubit. — The ground state that contains the re- 
sult of a given standard algorithm is specified as follows 
[l^ . First consider a particularly simple computation 
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FIG. 1: Electron in array of quantum dots performing a sin- 
gle qubit, 4 step Deutsch-Jozsa algorithm for a function f{x) 
where /(O) = /(!)• Qubit starts in logical 0, undergoes a 
Hadamard gate, the mapping |a::) (-l)^(^' |j;>, and a final 
Ifadamard gate. Each of the 4 columns of dots corresponds to 
an algorithm step. Near (far) dot of each column corresponds 
to logical (1). a) When AQC starts, Ifamiltonian parameter 
A = 0. Ground state is localized on logical dot of input col- 
umn, b) As we increase A, potential of first columns increases 
and tunneling between columns turns on. Ground state spills 
from input step into later steps, c) When A = 1, ground state 
amplitudes in final column carry calculation results. 




involving only a single qubit with basis states |0} and 
In the circuit model the qubit evolves through N + 1 
time steps: its initial state and a state after each algo- 
rithm step. If the initial state is |0) and algorithm step i 
consists of application of a 2 x 2 unitary gate Ui, then the 
two amplitudes at time step i are given by Ui ■ ■ - UilO), 
where 1 < i < A^. Since there are 2 amplitudes at each of 
the iV-l- 1 steps, the whole trajectory can be described by 
giving 2(A^-|-1) complex amplitudes. In GSQC, instead of 
a time-dependent state in a 2-dimensional Hilbert space, 
the qubit has a time- independent state in a 2(A^ + 1)- 
dimensional Hilbert space, with basis states c|q |vac) and 

c|]^|vac), i = 0,...,N. Here, {ci,x) is a fermionic 
creation (annihilation) operator for a particle in state 
X G {0,1} in mode i G {0, l,...Af} [l^. The ampli- 
tude of the time-independent wavefunction in basis state 
c|q I vac) {cl ^ I vac)) contains the amplitude of the time- 
dependent system in the basis state |0) (|1)) after algo- 
rithm step i. 

Illustrating this with a concrete physical system, 
Fig. Ic shows a 2-dimensional array of quantum dots (i.e., 
potential minima) with 4 columns of 2 dots each. Column 
i contains a state localized on the near dot (c| q |vac)) and 

a state localized on the far dot (c| ^ |vac)); this sample al- 
gorithm has -I- 1 = 4 steps. 



It is convenient to group creation operators into row 
vectors C| = cl o "^I i ■ Then, the (unnormalized) 
ground state containing the results of the algorithm is 
1*"^) = (^0 [U+Clu, [1]+---+cI,Un---U, [J])|vac). 
The results, stored in the states cjy q |vac) and cjy ^ jvac), 
can be extracted reliably 15(. To execute the GSQC, one 
realizes a specific Hamiltonian H{1) whose ground state 
is the time-independent state we have identified. When 
the computation involves a single qubit, the Hamiltonian 



takes the form H{1) = J2i=i h^{Ui) where 



(1) 



and where £ sets the energy scale. Here we have used the 
quadratic form CjVCj = J2x,y&{o,i}^^y4.x'^j,y^ where 
Vxy are the matrix elements of V. In the quantum dots 
example of Fig. Ic, H{1) controls the onsite energy of 
each dot and the tunneling coupling between each dot 
in one column and each dot in the next column. One 
can confirm |^^^ is an eigenstate with eigenvalue 
{H{1) IvP^) = 0); since every term ([1]) is clearly posi- 
tive semi-definite, IvP^) is the ground state |18|. 

We have a time- independent state that contains the re- 
sult of any given standard algorithm for one qubit and is 
the ground state of a known Hamiltonian H{1). Now we 
show that placing the system in this ground state can be 
done efficiently via the method of AQC, which constitutes 
a proof that AQC can simulate the circuit model (for a 
single qubit). To do so we introduce the Hamiltonian 
H{s) = J2Zi h'{X{s)U^) where A : s G [0, 1] [0, 1] such 
that A(0) = and A(l) = 1. If A = 0, then h'{0) = SCjC, 
reduces to a simple onsite energy term. There is then no 
tunneling from algorithm step i to algorithm step i + l. If 
A = 1, we recover H{1) = J2iLi h'^iUi) with h^{Ui) being 
the full operator llj. 

The (unnormalized) ground state of H(s) is simply 
I^E'^) as written above, but with a factor of A(s) in front 
of each unitary operator Ui. Alternatively, the ground 
state is given by the recursion relation 

1*^(5)) = (1 + Cj(A(s)t/,)Q_i) \^^-\s)) . (2) 

Intuitively, the state of a j step calculation, |5'^(s)), is 
formed by adding to the state of a j — 1 step calculation 
|\['^~^(s)) a term which annihilates the particle at j — 1 
and creates a particle at j with X{s)Uj applied to its state. 
The initial input state is simply |vl''^(s)) = Cq [^] |vac). 

Fig. la shows how the wavefunction is localized on the 
input row when A = 0. As expected from the form of 
h'^{\{s)Ui), as A increases, the onsite energy rises on rows 
0, . . . , A'^— 1, as in Fig. lb. The tunneling matrix elements 
in ll]) also begin to turn on, and the ground state wave- 
function starts to spill into the full array of dots. When 
A reaches 1, the wavefunction reaches the final row, as in 
Fig. Ic. 
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The traditional statement of the adiabatic theorem 
[lit is that the increase in A(s) must be sufficiently 
gradual that the system does not transition to an 
excited state as s goes from to f. If the time to 
take s from to 1 is T, transitions are suppressed if 
T ^ h maxs 



dHjs) 
ds 



where {s)} is any excited eigenstate of H{s), and 
Ej^{s) is its energy. Recent work has emphasized that 
this is not necessarily the right condition, since what 
really matters is not suppression of transitions through- 
out the entire adiabatic quantum algorithm, but rather 
that the overlap between the ground state of H{1) and 
the final adiabatic wavefunction be large (2^, |2l[. An 
adiabatic condition arises of the form T = hO{/S.E^^), 
where Ai?min is the minimum energy gap between the 
ground and first excited state, as s goes from to 1 
|20l . [2l| . Here we use this latter condition to prove that 
AQC can simulate the circuit model efficiently; we omit 
the (similar but more complicated) proof that uses the 
traditional adiabatic condition. We note, however, that 
we will need to use our knowledge of the gap structure 
(e.g., position of AE'min as a function of s) in order 
to achieve the running time T = hO{l^E^2^^)^ since it 
requires slowing the evolution down near AE'niin. 

To compute the minimum gap for an + 1 step cal- 
culation, we look for solutions of D'j^j^^{E) = 0, where 
Dl;^^ = det {H{s)/£ - E), H{s) is the 2{N + 1) x 
2{N + 1) Hamiltonian matrix, and E = E/£. We first 
make a unitary transformation to new operators Ci = 
{U}--- Ul)C^ which transforms H{s) = h'{X{s)U^) 
to X)t=i h^{^), where T is the 2x2 identity matrix. 
Writing out the matrix H{s), we find the iterative rela- 
tion Dn+i = (1 + - E)Dn - X'^Dn-1- The solutions 
to Dat+i — identify the exact single qubit eigenener- 
gies, which we find to be: Eq^s = and En^s — ^((1 — 
A(s))2 + 2A(s)(l - cos ^)) for n = I, . . . , A. By min- 
imizing the first excited state energy Ei_s with respect 



to A, one sees that Ei^s > iSsin 



(AT+l) 



£ 0(1/A2 



The minimum occurs when A = cos ■ Thus AE'min = 
Ei^s — Eq^s — £ 0(1/A^), and the simulation time is 
T = hO{AE^lj = {h/£)0{N^). Since this is polyno- 
mial we see that AQC can efficiently simulate the circuit 
model. 

Multiple qubits. — The recursion relation 
^ generalizes immediately to algorithms in- 
volving AI noninteracting qubits: |^'^(s)^ = 

n^=i(l + cI^{XUa^j)Ca^j-i)\^'-Hs)) where 
l^'O(s)) = n^tiC]^ [I] |vac). The multiple qubit 
Hamiltonian is just the sum of the single qubit Hamil- 
tonians H{s) = J2A=i^iLi ^AiH^WA.i)] one can 
verify that H{s) 1^-^(5)) = for arbitrary A. The AQC 
procedure for non-interacting qubits simply consists 
of the single qubit procedure applied to each qubit 
independently. 



Now, we allow the qubits to interact via two-qubit 
gates such as a controUed-NOT (CNOT). Suppose the 
algorithm specifies a CNOT gate between qubits A and 
B at line j. Then, instead of applying the factors 
(I + Cl^iXUA^,)CA^,-i) and (J + ^ /AC/b,,)Cb,,-i) 
to |*J-i(s)), we write 

l^'is)) = iI + c\^^,XcA^,,^,oCl^{XI)CB,j-i 
+ci,,,iAcA,,-i,iCL,,(Aa.)CB,,-i) \^'-'{s)) ■ 

If qubit A is in state 0, this operator applies (/ + 
Cg j{XI)CB.j-i) to qubit B. This is just the usual 
recursion relation factor that subjects qubit B to an 
IDENTITY gate. The factor for a NOT gate, (/ + 
Cg j{Xax)CB,j-i)), is applied to B if ^ is in state 1. 

When a CNOT gate is present, H{s) needs to be 
changed so that we still have H{s) |^'^(s)) = 0. One 
replaces terms h-'j^{XUA,j) and h-'g{XUB,j) in H{s) with 



h^iX, CNOT) = + hl,JN) 



(3) 



Here /i^b(ID) 



£{CB.jCA,j,0 

X^Cf 



X^Cbj-icaj-i.q) 



and h^A.sC^) = 



X {CB,jCA,jfl — X C_Bj_iCA,j-l,o), 
^{CB,jCA,j,l — X^axCB,j-lCA,j~l,l)'^ 

X {CB,jCA,j,i — X^axCB,j-iCA,j-i,i) are two-particle ana- 
logues of the one-particle IDENTITY gate h\{XJ) 
and NOT gate /i^Act^) defined by (P). The 



third term ^^^(P) 



^ J2i<j,k>j ^A. 



iCA,iCg^l.C, 



B,k 



C\ k'~'A,kCg iCB,i imposes an energy penalty on states 
in which one qubit has gone through the CNOT gate 
without the other. Since /i^ ^(A, CNOT) is positive semi- 
definite along with all the other terms in H(s), and as 
one verifies H{s) |^'^(s)) = 0, 1*^(3)) is stih the ground 



18| 



state 

To determine the effect of a CNOT gate on the gap, 
consider first a simple calculation with M = 2 qubits 
and a single CNOT at row j. Divide the Hamilto- 
nian into H = Hq + Hi, where Hi is just the CNOT 
gate ^ and Hq contains all of the single qubit terms. 
We know all of the exact eigenstates and eigenvalues of 
Hq from the single qubit analysis of Dat+i above. If 
the interaction Hi were absent, the qubits would sim- 
ply occupy these eigenstates of Hq independently. Let 
\Z) be a two-qubit state satisfying Hq \Z) = 0, where 
both qubits are in zero-energy ground states. Let 
be a state in which at least one qubit is excited. Our 
single qubit analysis of -Dat+i yields the exact result 
{Z\ Ho \Z) > £sm^ = £0{1/N^). 

The CNOT Hamiltonian Hi couples these eigenstates 
of Hq. It can be diagonalized analytically in the small 
basis of ground states \Z); for all states \Z) orthogo- 
nal to the computationally meaningful |^^(s)), we find 
{Z\H\Z) = {Z\Hi\Z) = £0{l/N'^). This is a rigor- 
ous variational upper hound of the true gap (because the 
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small number of \Z) states are not a basis for the full 
[2(A^ + l)]^-dimensional Hilbert space), and it is also a 
good estimate of the exact energy gap, as confirmed in 
the 2 qubit numerical calculation shown in Fig. 2. As a 
result, AQC requires a time T = {ni£)0{N'^). While this 
estimate is intuitively correct and numerically verified, 
what we need for our equivalence proof is a lower bound, 
which we now show is SO{l/N'^). To obtain this bound, 
note that an arbitrary excited state oi H = Hq + Hi can 
be written as a\Z) + (3\Z^ for some unique normalized 
\Z) and \Z). Since Hi is positive semi-definite, we find 

{Z\Hi\Z) {Z\Hi\Z) > \{Z\Hi\Z)\^. Given this inequal- 
ity, we find that the minimum of (H) with respect to a 
and f3 satisfies 

(H) > {Z\ H \Z) {Z\ Ho \Z) /{{Z\ H \Z) + {Z\ H \Z)). 

(4) 

We can estimate the numerator using the EO{l/N'^) 
bounds on (^Z\ Hq |Z) and {Z\ H \Z) stated above. Since 
the denominator of (|4]) certainly does not increase with 
N, the energy (H) is at least £0{l/N'^). 

A similar argument works even when the system has 
many qubits and many CNOT gates. We write the 
many-qubit Hamiltonian as H ~ Hq + Hi, where Hi 
includes all of the CNOT gates ^ and Hq contains 
all of the single qubit terms. The exact eigenstates and 
eigenenergies of Hq are immediately known from the sin- 
gle qubit analysis. An arbitrary state can be written 
a\Z) + P \Z'^ where in \Z) all qubits are in ground states 
oi Hq while in |Z) there is at least one excited qubit. We 
still have {Z\Ho\Z) > Ssin^jj^ = £0{l/N^). We 
can also show that {Z\ H \Z) > £0{1/N^) since for each 
term in \Z) there is always at least one CNOT gate that 
contributes to its energy. Using (|3]), we find the same 
£0{1/N'^) bound on (H). To facilitate extraction of the 
results of the computation, it is important that when the 
system is measured every qubit have a large amplitude on 
the final row N. As Ref. [IJ] shows, it is straightforward 
to adjust the Hamiltonian to ensure this happens. How- 
ever, the reduction of qubit amplitude at earlier stages of 
the calculation leads to a reduction of the gap. The esti- 
mate/upper bound becomes £0(1/N'^M) and the lower 
bound is £0{l/N^M'^). 

We have presented an explicit adiabatic procedure that 
will carry a system adiabatically into a ground state con- 
taining the result of an arbitrary standard quantum com- 
putation. Ref. [ll] shows how to use quantum telepor- 
tation to trade an arbitrary GSQC with steps and M 
qubits for a different GSQC with 7 steps, (2A^ - l)M 
qubits. Once the Hamiltonian is adjusted to facilitate 
extraction of the results [13], the running time of the 
new calculation is T = {Ti/ £)0{N'^AP). The upper 
bound/estimate of the gap yields T — [h/£)0{NM), 
which is proportional to the "volume" of the algorithm. 

We thank M.L. Cohen for useful discussions. We 
gratefully acknowledge support from the David and Lu- 
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FIG. 2: Minimum gap in units of £ for a 2 qubit system of 
TV steps. Results are shown for computations that entangle 
the two qubits into a Bell state and then disentangle then. 
A string of single qubit gates is also included at one of three 
stages of the computation. The minimum gap is roughly a 
linear function of l/N"^ . 
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